arXiv: 1508.0403Ivl [gr-qc] 17 Aug 2015 


Self-force as a cosmic censor in the Kerr overspinning problem 


Marta Colleoni, Leor Barack, Abhay G. Shah, and Maarten van de Meent 
Mathematical Sciences, University of Southampton, Southampton SOU IBJ, United Kingdom 

(Dated: August 18, 2015) 

It is known that a near-extremal Kerr black hole can be spun up beyond its extremal limit by 
capturing a test particle. Here we show that overspinning is always averted once back-reaction from 
the particle’s own gravity is properly taken into account. We focus on nonspinning, uncharged, 
massive particles thrown in along the equatorial plane, and work in the first-order self-force approx¬ 
imation (i.e., we include all relevant corrections to the particle’s acceleration through linear order in 
the ratio, assumed small, between the particle’s energy and the black hole’s mass). Our calculation 
is a numerical implementation of a recent analysis by two of us [Phys. Rev. D 91 , 104024 (2015)], 
in which a necessary and sufficient “censorship” condition was formulated for the capture scenario, 
involving certain self-force quantities calculated on the one-parameter family of unstable circular 
geodesics in the extremal limit. The self-force information accounts both for radiative losses and 
for the finite-mass correction to the critical value of the impact parameter. Here we obtain the 
required self-force data, and present strong evidence to suggest that captured particles never drive 
the black hole beyond its extremal limit. We show, however, that, within our first-order self-force 
approximation, it is possible to reach the extremal limit with a suitable choice of initial orbital 
parameters. To rule out such a possibility would require (currently unavailable) information about 
higher-order self-force corrections. 


I. INTRODUCTION 

Consider a Kerr black hole of mass M and dimension¬ 
less spin a/M = 1 — e^, with e ^ I; and send a particle of 
rest mass pL M into the black hole along the equato¬ 
rial plane. Ignore back-reaction effects on the particle’s 
motion, and assume it follows a geodesic of the Kerr ge¬ 
ometry, with conserved specific energy E and angular 
momentum L. Jacobson and Sotiriou [T] showed that 
there exists an open domain in the {pi,E,L'\ space, for 
which the final configuration is overextremal: 

aM + fj.L > {M + nEf. (1) 

Later work by Barausse et al. nig demonstrated that, at 
least for some orbits in that open domain, radiative losses 
cannot prevent the overspinning. The purpose of our 
current work is to show that the overspinning scenario is 
entirely ruled out once all relevant back-reaction effects 
are taken into account. 

The history of the problem goes a while back. In 1974 
Wald d] proposed the black hole-particle system as a 
test bed for weak cosmic censorship. He showed (ignoring 
back-reaction) that an initially extremal Kerr-Newman 
black hole cannot be made overextremal by capturing a 
particle, although parameters can be chosen so that the 
black hole remains extremal. Particles carrying “dan¬ 
gerous” amounts of angular momentum, electric charge 
and/or spin are never captured; rather, they get deflected 
away by centrifugal, electrostatic and/or spin-spin cou¬ 
pling forces. In 1999 Hubeny [S] argued that the situa¬ 
tion is different if one starts with a nearly extremal black 
hole: focusing on the problem of an electric charge in 
a Reissner-Nordstrom geometry, she showed there is an 
(open) set of captured orbits that overcharge the black 
hole in that case. However, she also reasoned that back- 
reaction from the electromagnetic self-force may well pre¬ 


vent suitable particles from being captured. 

Recent work has gone some way towards confirm¬ 
ing that expectation. Isoyama, Sago and Tanaka 
[B] considered an electric charge released from the 
point of quasi-static equilibrium—an initial configura¬ 
tion which they represented by the exact (static) “dou¬ 
ble Reissner-Nordstrom” solution, whose total Arnowitt- 
Deser-Misner (ADM) energy and charge are known an¬ 
alytically. Showing that radiative losses from the subse¬ 
quent plunge into the black hole are negligible, Isoyama 
et al. calculated that the final configuration cannot be 
that of an overcharged geometry. They further argued 
that the same conclusion must then apply to any radi¬ 
ally falling charge. In a later work, Zimmerman, Vega 
and Poisson [7] presented an explicit (numerical) calcu¬ 
lation of the charged particle’s trajectory including the 
full effect of the electromagnetic self-force. Analyzing a 
large sample of orbits within the domain identified by 
Hubeny, they found no example of successful overcharg¬ 
ing: All particles with charge and energy suitable for 
overcharging the black hole were found to be repelled 
before reaching the horizon. 

However, both above analyses 0 [7] have neglected 
the coupling to gravity (both the back-reaction from the 
gravitational perturbation sourced by the perturbed elec¬ 
tromagnetic stress-energy, and the perturbation in the 
electromagnetic field due to the particle’s mass), which 
cannot be easily justified for the problem at hand. A 
complete analysis would require a calculation of both 
electromagnetic and gravitational self-forces in the cou¬ 
pled problem, which remains a difficult challenge despite 
recent progress [HHin]. 

Here we instead consider the purely gravitational prob¬ 
lem of an electrically-neutral particle in Kerr geometry. 
Thus evading the coupling problem comes at the obvious 
cost of giving up the convenience of working in a spheri- 
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cally symmetric background. However, the current state 
of affairs in self-force calculations is such that computa¬ 
tions on a Kerr background are now routine (see, e.g., 
the recent my Here we revisit the problem from this 
new vantage point, and provide what we believe to be 
a first example of a complete, self-consistent analysis of 
the overspinning problem within the first-order self-force 
approximation. 

The necessary groundwork for our calculation was laid 
in a previous paper by two of us [T^], to be referred to 
in what follows as CB. Focusing on equatorial capture 
trajectories, and excluding deeply bound orbits (see be¬ 
low), CB formulated a set of “censorship conditions” that 
must be satisfied in order for the overspinning scenario 
to be ruled out for all such orbits. The conditions are 
formulated in terms of a certain one-parameter family of 
geodesic trajectories, and involve the gravitational self¬ 
force (GSF) evaluated along the trajectories in the limit 
of an extremal geometry. In the current work we numer¬ 
ically evaluate the censorship conditions, and establish, 
with confidence, that they are indeed satisfied. 

Our results reveal some interesting nuances. To de¬ 
scribe them we must first review the results of CB in 
some more detail. CB first identify (working at leading 
order in e ^ 1 ) the complete region in the parameter- 
space {E,L,e,r] := fi/M} for which Eq. Q is satisfied, 
i.e. overspinning occurs, if the GSF is ignored and the 
particle follows a geodesic of the Kerr background. It is 
found that, for any E > 1, overspinning occurs in cer¬ 
tain narrow ranges (oc e) of 77 , L values. Geodesics with 
E < 1 cannot overspin. Focusing, therefore, on particles 
sent in from infinity, CB proceed to incorporate all rele¬ 
vant GSF effects through 0{rf) in Eq. Q. This includes 
gravitational-radiation losses, as well as the (more subtle 
but not less important) 0{r]) correction to the critical 
value of impact parameter for capture. Since the cen¬ 
sorship condition is formulated in terms of background 
(“black-hole centred”) coordinates, which is beneficial in 
practice for GSF calculations, rather than in a “center-of- 
mass” system where the identification of ADM quantities 
is easier, care is taken in expressing relevant ADM quan¬ 
tities in terms of background-related quantities through 
the required order [0{rf‘) in Eq. §] . CB thus obtain an 
inequality involving the initial energy and angular mo¬ 
mentum (at infinity), as well as 77 , e and the GSF, which 
describes a condition for the final configuration to remain 
subextremal. 

CB then show that the above censorship condition 
is satisfied for all capture trajectories (and for all suf¬ 
ficiently small 77 , e) if and only if it is satisfied for the 
one-parameter family of “critical” orbits—those lying on 
the scatter-capture separatrix in the space of initial con¬ 
ditions. Thus, the problem may be conveniently reduced 
to the question of whether critical orbits satisfy the cen¬ 
sorship condition. CB parametrize such orbits by the 
energy-at-infinity, and formulate a necessary and suffi¬ 
cient censorship condition on that one-parameter family. 

The behavior of critical orbits is subtle when radiation 


is taken into account. CB distinguish between “strong” 
(exponential) and “weak” (power-law) fine-tuning of the 
initial conditions, depending on the proximity to exact 
criticalityQ Strongly fine-tuned orbits radiate a signifi¬ 
cant fraction of their initial energy as they evolve adia- 
batically through a sequence of unstable circular orbits; 
in the most extreme case (“perfect” fine-tuning) “all” the 
orbital energy is radiated away and the particle plunges 
upon reaching the innermost stable circular orbit (ISCO) 
[13] . Weak fine-tuning still results in trapping on quasi¬ 
circular orbits, formally for an infinite amount of time as 
77 —>■ 0 , but the fraction of energy emitted is vanishingly 
small in that limit. (A more precise definition of weak 
and strong fine-tuning will be given in Sec. HH) 

CB show that, in the case of weak fine-tuning, all 
radiation-related terms drop out of the censorship condi¬ 
tion through the relevant order [ 0 ( 77 ^)]. Only one GSF 
term remains, associated with the conservative correction 
to the critical value of the angular momentum (at fixed 
initial energy). CB describe two methods for calculating 
this correction. The first, direct method involves integra¬ 
tion of the GSF along the critical geodesics all the way in 
from infinity and down to the limiting (unstable) circular 
orbit. The second method, which relies on a recent per¬ 
turbative formulation of the conservative Hamiltonian for 
circular orbits [HHisi, requires information about the lo¬ 
cal metric perturbation (not the GSF itself) only on the 
limiting circular orbit. The second method is much more 
easily implemented, and will be used here to compute the 
shift in the critical value of the angular momentum. 

In the case of strong fine-tuning, radiative terms enter 
the censorship condition already at the leading order in 
77 . GB show that the additional information required in 
this case amounts to a single function of the circular- 
orbit energy, namely the ratio Tl{E) between the flux of 
energy absorbed by the black hole and the flux of energy 
radiated to null infinity. 

Thus, two separate censorship conditions emerge. The 
first applies to weakly fine-tuned orbits and involves only 
information about the conservative piece of the GSF. It 
has a remarkably simple form, especially in conjunction 
with the Hamiltonian formulation of [HHis]. The sec¬ 
ond censorship condition applies to strongly fine-tuned 
orbits, and requires also knowledge of the radiative flux 
ratio TZ{E). Both conditions involve quantities that are 
formally evaluated in the limits e —>■ 0 and 77 —> 0. Thus, e 
and 77 themselves do not features in the final form of these 
conditions. The censorship inequalities should be inter¬ 
preted as necessary and sufficient conditions for overspin¬ 
ning to be averted for sufficiently small e and rj. However, 
there is no assumption about the relative magnitudes of 
these two parameters. 


^ In CB’s terminology, these two types are respectively referred to 
as “fine-tuned” and “generic” near-critical orbits. We prefer here 
“strong” and “weak” fine-tuning, to avoid possible confusion. 
Both types are fine-tuned, but at a different level. 
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In the current paper we evaluate both censorship con¬ 
ditions. We consider weakly fine-tuned orbits first. We 
give a strong numerical evidence to suggest that all such 
orbits precisely saturate the overspinning condition, i.e. 
they lead to a precisely extremal final geometry, for any 
value of the initial energy. Next we examine whether 
over-extremality may be achieved through a strong fine- 
tuning of the initial conditions. This is an intriguing pos¬ 
sibility, not least because violation of censorship has been 
(famously) demonstrated in the past for other carefully 
fine-tuned configurations m- However, we find that 
strong fine-tuning in our problem always acts to drive 
the system away from extremality. Capture of strongly 
fine-tuned orbits results in a subextremal geometry. 

A corollary of the above results is that, within our first- 
order GSF approximation, captured particles generically 
leave the black hole subextremal, except in the case of 
weakly fine-tuned critical orbits, which appear to pre¬ 
cisely saturate the censorship condition. In the latter 
case, one must go beyond the first-order GSF approxi¬ 
mation to determine whether saturation actually occurs 
in the physical problem, and theoretically there even re¬ 
mains the possibility of overspinning. One of our main 
conclusions is that a definitive answer to the overspinning 
question cannot be reached within the first-order GSF ap¬ 
proximation, just as it could not be reached within the 
geodesic approximation. This is perhaps disappointing, 
but also intriguing. 

The rest of the paper is structured as follows. In Sec- 
tion|TT]we give a more detailed review of the GB analysis, 
and present the censorship conditions. In Section [lll| we 
evaluate the condition in the case of weak fine-tuning 
and conclude that it is saturated, but overspinning never 
occurs. In Section |IV] we show that strong fine-tuning 
promotes censorship in our problem. Section [V] contains 
a summary and a discussion of our results and their im¬ 
plications. 

Throughout this paper we set G = c = 1 and use the 
metric signature (—, -|-, -I-, -b). 

II. REVIEW OF CB: THE CENSORSHIP 
CONDITIONS 

A. Basic setup 

The system in question is a gravitationally-bound bi¬ 
nary of mass ratio rj := p./M 1 . The larger body is a 
near-extremal Kerr black hole of spin aM = M^(l — e^), 
with e <C The smaller body (“particle”) is a com¬ 
pact object of negligible spin, which, for concreteness and 
simplicity, we may think of as a small non-spinning black 
hole. We work in black-hole perturbation theory, which is 


^ Beware the alternative convention a/M = (1 — ~ 1 — e^/2 

is common in related literature. 


based on a formal expansion in 77 about the background 
geometry of the large black hole. Specifically, we will 
consider the binary dynamics in the first-order GSF ap¬ 
proximation, namely through the first order in rj beyond 
the limit of geodesic motion. We let {t,r, 9 ,(p} denote 
standard Boyer-Lindquist coordinates on the background 
geometry, and let the particle be sent in along the equa¬ 
torial plane, 9 = 7r/2. From symmetry it is clear that the 
orbit will remain equatorial even under the effect of the 
GSF. 

It is assumed that the standard GSF formulation of 
the binary dynamics applies: the particle’s motion is de¬ 
scribed locally in terms of an effective (accelerated) tra¬ 
jectory on the background spacetime, which, at our order 
of approximation, is insensitive to the particle’s internal 
structure. In particular, the particle will be assumed to 
have fallen into the large black hole if and only if its 
trajectory had crossed the latter’s event horizon. 

It will be useful to define, for given {M,/r, e}, 

E := (Madm - M)/p, (2) 

L := (Jadm - aM)/{pM), (3) 

where Madm and Jadm are the total (conserved) ADM 
mass and angular momentum of the binary’s spacetime. 
E may be interpreted as the particle’s contribution (per 
p.) to the ADM energy in the “black hole frame”, and L 
is the particle’s specific (and a-dimensionalized) orbital 
angular momentum, again in the “black hole frame” 
For the GSF formulation to make sense, we require, in 
addition to 77 ^ 1 , also E ^ I/77. We do not consider 
ultrarelativistic particles with E ^ 1 (not even ones with 
E <C l/v), because for such orbits there does not yet ex¬ 
ist a rigorous GSF formulation. Nor do we consider the 
family of “deeply bound” orbits discussed in Sec. II.B of 
GB: low-energy orbits that are confined to the immediate 
exterior of the large hole, below the ISGO radius. GB ex¬ 
plain that, for e ^ 1 and relevant values of 77, such orbits 
plunge into the black hole within an amount of proper 
time comparable to the particle’s own “light-crossing” 
time. Ignoring the particle’s finite extent does not seem 
justified in that situation, so such orbits require a sep¬ 
arate treatment. For simplicity, GB opted to exclude 
deeply bound orbits from the analysis. In practice, this 
amounts to assuming that the particle is thrown in from 
outside the ISGO radius. 

Our system is thus parametrized by the quartet of val¬ 
ues {77, e, E, L}, up to an initial orbital phase and up to 
a trivial overall scale set by M, neither of which are of 
relevance in our analysis. This parameter space splits 
into two precisely complementary subspaces, one corre¬ 
sponding to orbits that end up crossing to horizon, and 
another corresponding to orbits that scatter to infinity 
without crossing the horizon. For captured orbits, let 


® We use E and L in place of CB’s and I^adM’ rotational 

simplicity. 
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Mb and Jb be the Bondi energy and angular momen¬ 
tum of the spacetime at the (retarded) time of horizon 
crossing. The question of overspinning then takes the fol¬ 
lowing simple form: Among all configurations {77, e, E, L}, 
with arbitrarily small {77, e}, are there any in which the 
particle crosses the horizon and Jb > Mg? If so, then the 
likely scenario is that a naked singularity is exposed [1]. 
If not, then the conjecture of cosmic censorship holds. 

B. Test-particle limit 

Let us first review the situation when back-reaction is 
ignored and the particle is treated as a test mass (this is 
the case first analyzed in [T], but the details below follow 
CB). The particle moves on a timelike geodesic of the 
Kerr background metric with conserved (specific) 
energy and angular momentum given by E = —ut and 
L = u^jM, where u°‘ is the particle’s four-velocity and 
Ua = ■ Assuming the particle is thrown in from 

outside the ISCO radius, it will cross the horizon if and 
only if L is smaller than a certain critical value Lc°^(E; e), 
which corresponds to the angular momentum of an un¬ 
stable circular geodesic orbit with energy E (“peak of 
the effective potential”). We use a superscript ( 0 ) to de¬ 
note the test-particle limit. CB find, in the near-extremal 
case, 

L(°)(E; e) = 2 E + (hE^ - 2 Y/‘^e + ( 4 ) 

The Bondi mass and angular momentum are constant 
and given by Mb = M -t- /iE and Jb = aM + fiML, re¬ 
spectively. Captured orbits overspin the black hole if and 
only if they satisfy Jb > Mg, which may be rearranged 
to give L > 2E-|-e^/77-|-?7E^. Thus, for given {g,e, E}, the 
capture condition bounds L from above, L < while 
the overspinning condition bounds it from below. Ne¬ 
glecting high-order terms in e, overspinning would occur 
if the double inequality 

eV ?7 + ? 7 E^ < L- 2 E < ( 6 E 2 - 2 )i/ 2 e ( 5 ) 

can be satisfied with e, 77 ^ 1. 

A simple analysis [T^] shows that ([^ is never satisfied 
for E < 1 , but it can be satisfied for any E > 1 by choosing 
77 from within the range 

e/_(E)<77<e/+(E), (6) 

with f± = [\/ 3 E 2 — 1 ± VE^ — 1]. For any such 

choice, there is an open range of L values satisfying ([^. 
Thus, for a given e ^ 1 , there is an open overspinning 
domain in the space of {77, E, L}, described by E > 1 
together with Eqs. § and (§. This domain is depicted 
in Figures 2 and 3 of CB. The key point: if back-reaction 
effects were negligible, it would be possible to overspin 
the black hole by throwing a particle in from infinity. 

Equation Q shows that, as might be expected, 77 and 
e are of the same order of magnitude in overspinning 


configurations. Equation ([^ shows, in turn, that all 
overspinning orbits have L — 2 E = 0 {g). Such orbits 
are located close to the parameter-space separatrix be¬ 
tween captured and scattered orbits. They each exer¬ 
cise 0(ln77) 3 > 1 near-circular revolutions near the peak 
of the effective potential before plunging into the black 
hole. The amount of excess spin produced by such orbits 
is Jb - M| = 0 {rf). 

This implies that back-reaction effects may qualita¬ 
tively change the outcome of the analysis. A priori, 0(77) 
self-acceleration may shift the value of the critical angular 
momentum (for a given energy) by an amount of 0(?7), 
comparable to the width of the overspinning window. 
Furthermore, radiation of energy and angular momen¬ 
tum in gravitational waves to infinity may, a priori, alter 
the final excess spin by an amount of 0(77^ In 77), greater 
than the excess spin predicted with radiation neglected. 
Clearly, the point-particle approximation is inadequate in 
our problem. It is a case where a leading-order perturba¬ 
tive treatment predicts its own inadequateness. By the 
end of our analysis it will become clear that, remarkably, 
the same conclusion carries over to the next perturbative 
order. 


C. Self-force approximation and the correction to 
the capture—scatter separatrix 

In the first-order GSF approximation, the particle’s 
motion is described locally by an accelerated world¬ 
line in the Kerr background metric gap- The tangent 
four-velocity, normalized using gapu°‘u^ = —1, satisfies 
/i7t^V/37t“ = F“, where is a covariant derivative com¬ 
patible with gap, and F“(oc 77^) is the GSF. The latter 
may be attributed to a certain smooth (locally-defined) 
self-potential denoted —the “R field”—which is a 

particular solution of the source-free linearized Einstein’s 
equations [ 18 ]. The particle’s trajectory may also be in¬ 
terpreted as a geodesic in gap + h^p. 

For given e, 77, our equatorial orbits may again be 
parametrized by the ADM-related constants E, L defined 
in Eq. (§. These now also carry information about the 
self-gravity of the particle and about the radiation con¬ 
tent of spacetime. For given e, 77, E there exists a crit¬ 
ical value Lc(E;e,77) representing a threshold for imme¬ 
diate capture: orbits with L > scatter off the black 
hole (at least at first approach; they may still eventually 
fall into the black hole following subsequent encounters), 
while those with L < L,. are immediately captured. The 
geodesic limit, 77 —)• 0 , of Lc(E;e, 77) is Li°^(E;e), given in 
Eq. @. Not surprisingly, much of the physics relevant 
to our problem is played out near the capture-scatter 
threshold. 

In the geodesic case (and for a fixed e), L = Li°^(E; e) 
describes a 1-parameter family of disjoint critical orbits. 
Each critical orbit is homoclinic in nature, approaching 
an (unstable) circular geodesic at t —)■ 00. The situation 
becomes more subtle when GSF effects are included, due 
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to dissipation. For given e,r], the function L = Lc(E;e) 
still describes a 1-parameter family of critical orbits, but 
these are no longer disjoint. Rather, they all merge to 
form a “global attractor”. The attractor may be thought 
of as a smooth sequence of quasi-circular unstable or¬ 
bits starting at the “light ring” and ending at the ISCO. 
A critical orbit meets the attractor at a point that de¬ 
pends on its initial energy, and proceeds to evolve ra- 
diatively along it until it transits to a plunge near the 
ISCO. A small perturbation away from the critical value 
L = Lc(E;e, 77) results in the particle’s leaving the at¬ 
tractor before the ISCO is reached. The point of depar¬ 
ture may be controlled by fine-tuning the magnitude of 
L around its critical value. 

When considering near-critical orbits in the overspin¬ 
ning analysis, it is important to distinguish between two 
cases, differing by how well L is tuned to its critical value. 
“Strong” fine-tuning is one in which L — Lc ^ exp(—a/ry), 
with some positive constant a. Such orbits radiate away 
amounts of 0 {r]) energy and angular momentum as they 
evolve along the attractor [leading to 0(1) changes in 
the values of the particle’s specific energy and angu¬ 
lar momentum]. “Weak” fine-tuning is one in which 
L — Lc ~ T]^ , with some /? > 1 . Such orbits radiate 
away only 0(77^ In 77) of energy and angular momentum 
before leaving the attractor. As will be discussed below, 
CB found that radiation effects enter the overspinning 
condition at relevant order only for strongly fine-tuned 
near-critical orbits. For all other orbits, radiation terms 
are negligible within the first-order GSF approximation. 

The overspinning analysis requires an explicit expres¬ 
sion for Lc(E;e,77) through 0 {r],e). For orbits that come 
from infinity (i.e., ones with E > 1 ) CB obtain, to the 
required order, 

U(E;e, 77 ) = L( 0 )(E;e) + 5 U(E; 77 ), ( 7 ) 

where is given in Eq. Q, and the GSF correction is 

5 L, = -77(E2 + 1) + hm / ( 2 Ft + drjuF (8) 
<=^0 \ / 

Here Ft := Ft /and F^ := F^/{Mfi) are components of 
the self-acceleration (oc 77), and the integration is carried 
out along the critical geodesic with energy E and angular 
momentum Lc°^(E;e). The cutoff radius r = i?j(E;e) is 
that of the corresponding unstable circular geodesic. The 
term —77(E^ -|-1) arises from the transformation between 
the background-defined energy-at-infinity and angular- 
momentum-at-infinity, —Ut(r —)■ 00) and u^(r —>■ 00) re¬ 
spectively (these are used in the GSF description of the 
motion), and the ADM-related quantities E and L—see 
Eq. ( 72 ) and Appendix A of CB. 


D. Censorship conditions 

Consider an orbit that ends up crossing the event hori¬ 
zon, and let £+ and £+ be the total energy and angular 


momentum in gravitational waves radiated out to null 
infinity up until the retarded time of horizon crossing. 
The Bondi mass and angular momentum correspond¬ 
ing to that retarded time are Jb = M + — and 

aM -I- /rL — £+, respectively. The condition for overspin¬ 
ning is Jb > M^, which, upon substituting a/M = l — e^ 
and rearranging, reads 

-h 77(2E - L) - W+ + (77E - £:+)2 < 0 , ( 9 ) 

where W’*' := 2 £’+ — £+. CB showed that the left- 
hand side of this inequality is minimized (with respect 
to L, at fixed E, e, 77) by near-critical orbits with L = 
Lc(E; e, 77) -|- 0(77). Thus it suffices to restrict attention to 
this type of orbits. If, for some E (and e,r] 1 ), the 
inequality ([^ can be shown to apply for a near-critical 
orbit, that would establish a violation of censorship. If, 
conversely, can be shown not to be satisfied even for 
near-critical orbits (for any E), then it cannot be satis¬ 
fied for any captured orbit, and censorship is protected. 
We can therefore proceed by substituting L = Lc(E;e,77) 
in Eq. (©• Using 0 and 0 one thus arrives at the 
censorship condition 

-f er,(t){E) - pSldE) - W+ -b (77E - £+f > 0 , ( 10 ) 


where 


^(E) :=-(6E2-2)1/2. 


( 11 ) 


Censorship is protected if and only if ( | 10 [ ) is satisfied for 
all E when e, 77 are arbitrarily small. 

CB next showed that ( 10 ) can be written in the con¬ 
venient equivalent form 

e2 + e770(E)-77<5L™-(E)-W(+^) + (77E-£:+^))2 > 0, (12) 

where is the contribution to 5 \-c from the conser¬ 
vative piece of the GSF, and are the contribu¬ 

tions to ,£~^ from the quasi-circular part of the near- 
critical orbit. [CB established this by showing that the 
dissipative contribution to ^Lc cancels against the part of 
W*" corresponding to the “approach” part of the critical 
orbit; that the contribution to £'^ from the approach is 


negligible in Eq. (lOl; and that the contributions to both 


yV"*" and £■+ from the final plunge from the attractor to 
the horizon are also negligible in that equation.] In Eq. 
(12), is an 0(77) quantity for any near-critical or¬ 
bit, while the 77-scaling of the radiative terms and 

^(qc) depends on the degree of fine-tuning. For weakly 
fine-tuned orbits CB find = 0(6)0(772 In 77) and 

^(qc) “ O In p), while for strong fine-tuning the scal¬ 
ing is instead = 0 {e) 0 {r]) and = 0(77). 

It follows that, for weak fine-tuning, both radi ative 
terms and subdominant in Eq. ( 12 ) and 

drop out at our working order (we assume here e]m77| ^ 
1 ). The two terms survive in Eq. ( 12 ) only for strongly 
fine-tuned orbits. We proceed by considering the two 
cases separately. 
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1 . Weak fine-tuning 


In the case of weak fine-tuning we are left with the 
simple censorship condition 

-I-? 7 e(/)(E)-I-? 7 ^V'(E) > 0, (13) 

with ip{E) := — (5L(;°"®(E). We have introduced 




:= T] di; 


(14) 


so that both coefficients (^(E) and '0(E) in Eq. (13) have 
finite (generally nonzero) limits rj ^ 0 and e —>■ 0 (taken 
with a fixed E). Censorship holds if and only if, for any 
E, the inequality in Eq. (13) is satisfied for all (small, 
positive) 77 , e. A simple analysis shows [12] that this de¬ 
mands 0 > 0^/4, leading to the necessary and sufficient 
censorship condition 


^Lco„s(E)< 1 (i_e2), 


(15) 


for all E. The condition is sufficient in the sense that 
its confirmation would establish that overspinning is not 
possible for any rj,e 1 and any E (when strongly fine- 
tuned orbits are excluded). The condition (15) is nec¬ 


essary in the sense that its violation for some E would 
imply that, for that value of E, there exist 77 , e ^ 1 with 
which overspinning can be achieved. 

In Eq. (15) we may allow E to vary in the full range 
Eisco < E «C 1 / 77 , where Ejsco = ^ the ISCO value 
of E in the limit 77 , e —>■ 0. This is notwithstanding the 
fact that the explicit expression given for dL™''® in Eq. 
^ (with the full GSF Fa replaced with its conservative 
piece) only applies for E > 1. In subsection HE below we 


will derive an alternative expression for dL™"®, applicable 
for any E in the above full range. 


2 . Strong fine-tuning 


First, however, let us formulate a censorship condition 
for strongly fine-tuned orbits. As mentioned, in that case 
one has = 0 (e) 0 ( 77 ) and = 0 ( 77 ) and both 

terms feature already at leading order in Eq. (12). We 
introduce the rescaled quantities 


W(tc) 


:= (776)- 




£ + 
■^(qc) 


-V ^(qc)> 


(16) 


which should have finite (generally nonzero) limits rj^e ^ 
0. We note that and depend not only on E 

but also on the precise fine-tuning. It is then convenient 
to re-parametrize the problem using the pair {Ei, Ef} := 
{E, E — f/ 77 }, where £ is the total energy radiated away 
(both to infinity and down the black hole) during the 
quasi-circular whirl. To leading order in 77 , Ei is the 
particle’s specific energy just upon entering the whirl, 
and Ef is its specific energy just before leaving it. The 
difference Ei — Ef is the total energy (per 77 ) radiated 


during the whirl. Any value oi Ef in the range Eisco < 
Ef < Ei may be obtained via strong fine-tuning. 

The censorship condition (12) now takes the form 


77e0(£'j, Ef) -\- ?7^0(Ei, Ef) > 0, 


(17) 


with 

0 =-( 6 ^ 2 - 2 ) 1/2 

0 = -dLr^(E;0 + (E,-f+c))"- (19) 

Here we have used the fact that E = Ei to leading order 
in 77 . 0 and ■0 depend on Ef through kV(|jc) '^(qc)’ 
spectively. Censorship holds if and only if, for any Ei, Ef 
satisfying Eisco < Ef < Ei, the inequality in Eq. 0 is 
satisfied for all (small, positive) 77 , e. A simple analysis 
shows m that this happens if and only if 


■0 > ^min{ 0 / 2 , 0 }^ 


( 20 ) 


for any Eisco < Ef < Ei. This constitutes a necessary 
and sufficient censorship condition for strongly fine-tuned 
orbits. 

The evaluation of (20) requires, in addition to dL™"®, 
also the radiative quantities ^^jc) show 

that, at the required order, these can be conveniently 
obtained using 


dE 


^(qc)('®»'®/) l + n{E)' 


>V0 dEi,Ef) = -[ dE, 

(qc)V f! 1 + TZ{E) 


where 


b{E) := 6E(6E2 - 2)"1/2. 


( 21 ) 

( 22 ) 

(23) 


Here Tl{E) is the ratio between the flux of energy ab¬ 
sorbed by the black hole and the flux of energy radiated 
to infinity, for a particle on a circular geodesic orbit with 
specific energy E, in the extremal limit e —>■ 0. Hence, the 
only information we require about radiation is encapsu¬ 
lated in a single dimensionless function TZ{E), evaluated 
on circular geodesics. We note TZ < 0 for E < ^, the 
superradiant regime (we will show how this result is ar¬ 
rived at in Sec. IV ). However, we also have TZ > —1 for 
any E, implied by the known non-existence of “floating” 
orbits. Since the integrand in (22) is positive definite 
(and Ei > Ef), we have VV(^c) ^ 9, so the sign of 0 in 
Eq. (20) is not a priori known. 


E. Reexpressing in terms of redshift 

Calculating dL™''® using Eq. has two drawbacks. 
First, this formula applies only to particles thrown in 
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from infinity (i.e., ones with E > 1). This restriction 
comes from the fact that in deriving ([^ CB relied on hav¬ 
ing at hand an explicit relation, correct through 0{r]^), 
between the ADM mass and angular momentum on one 
hand, and background-defined quantities like the four- 
velocity on the other hand. Such a relation could only 
be derived under the condition that the initial binary sep¬ 
aration was infinitely large. Although the test-particle 
analysis provides some motivation for concentrating on 
such orbits, it is preferable to relax the restriction E > 1 
when considering the GSF case. 

The second drawback of the formula Q is a practical 
one. Implementing the formula requires a calculation of 
the GSF along unbound orbits. However, existing GSF 
calculation methods and working codes are specialized 
to bound, quasi-periodic orbits. Adapting these codes to 
deal with unbound orbits might be possible in principle, 
but would require much development of new method and 
code. 

CB derived an alternative formula for circum¬ 

venting both problems. Their derivation built on recent 
work by Isoyama, Le Tiec and collaborators [HHm, in 
which expressions were derived through 0{rf') for the 
(Bondi-like) energy and angular momentum of a circu¬ 
lar binary with time-symmetric boundary conditions (so 
that spacetime admits a global helical symmetry). These 
expressions were given in terms of the so-called “redshift” 
variable [15] z := 1/u*, where is the four-velocity 
of the circular orbit, normalized not in the background 
metric ga^ but rather in the smooth perturbed metric 
gap + CB related Refs. [HHUl’s notions of energy 
and angular momentum to E and L, and thereby obtained 
an expression for in terms of the redshift z. 

The expression they derived is remarkably simple: 


5L" 


'(E) = -liin(5z(E;e), 

e—vO 


(24) 


where 5z is the first-order GSF correction to z, defined 
through the expansion 

z(E; e) = z(°)(E; e) + (5z(E; e) + 0{rf ). (25) 

Here E is used to parametrize the circular orbit, is 
the geodesic limit of z, and the 0 ( 77 ) GSF correction 
i5z is defined with a fixed E (and a fixed <!)0 We note 
that, at our order of approximation, it is permissible to 
evaluate Eq. (241 using a sequence of geodesic circular 
orbits, replacing the ADM-related quantity E with the 
geodesic specific energy E = —Ut- 

The relation (241 is valid for all E > Ejsco and can be 
used to derive for any unstable circular orbit, with¬ 

out the restriction E > 1. Furthermore, the evaluation of 
^Lcons requires only circular-orbit GSF data. 


^ Beware that, in related literate, the GSF correction to z is often 
defined with a fixed orbital frequency Q, not with a fixed E as 
here. 


readily computable using existing codes. In fact, using 
Eqs. (24), (21) and (22) it is evidently possible to evalu¬ 
ate both censorship condition (15) (for weak fine-tuning) 
and (20) (for strong fine-tuning) using circular-orbit in¬ 
formation only. Our calculation in the next two sections 
takes a full advantage of that fact. 


III. EVALUATION OF CENSORSHIP 
CONDITION 


In this section we evaluate the censorship condition 
( [T^ , which ignores the possibility of strong fine-tuning. 
The latter will be considered in Sec. El At our order of 
approximation, we may replace the ADM-related quan¬ 
tity E in Eq. (15) with the geodesic specific energy E, and 


regard (151 as a condition on the family of (unstable) cir¬ 
cular geodesic orbits, evaluated in the limit e —>■ 0. If the 
inequality can be shown to hold for all E > Eism = 
then overspinning is ruled out for all orbits (except, pos¬ 
sibly, strongly fine-tuned ones). 

The evaluation of (15) requires only the function 
jLcons(^) _ JL(;°"®(if)/? 7 , and we shall use Eq. (24) to 


calculate it. As will be shown, the evaluation of the 


redshift correction Sz in Eq. (24) becomes particularly 


simple in the limit e —S^ 0, to the effect that the essential 
part of the calculation can be done analytically. The only 
numerical input we shall require is a verification that a 
certain perturbative quantity has a finite limit e —>■ 0; the 
precise numerical value of that limit will not be impor¬ 
tant to us. Below we first present the analytical part of 
the calculation, and in subsection (HI BI we discuss the 
numerical input. 


A. Analytical considerations 


Isoyama et al. m show that the first-order GSE cor¬ 
rection to the redshift z can be obtained via 

Sz = -z‘'°'>H^ where := . (26) 

Here z^°^ (recall) is the geodesic limit of z (taken with 
fixed E,e), u°‘ is the four-velocity of the corresponding 
circular geodesic, and is the “time-symmetric” 

part of the R field evaluated on that geodesic. More pre¬ 
cisely, is a certain regular piece not of the physical, 

retarded metric perturbation, but of a time-symmetrized, 
“half-retarded plus half-advanced” perturbation that is 
responsible for the “conservative” part of the dynamics. 


The evaluation of JL™"® via Eq. (24) requires taking 
the limit e —>■ 0 of 5z with fixed E. To evaluate the limit 
of the factor z^^\ start with the general expression m 


= (1 - aD)^/2 [1 + _ 3(MD)2/3(I - aD)!/^ 


-, 1/2 


( 27 ) 




















in which := d(j)/dt = + is the orbital 

angular velocity, with R being the Boyer-Linquist radius 
of the orbit. The latter admits the small-e, fixed-i? ex¬ 
pansion [T^ 

R/M = l + epi{E) + e^P2{E) + 0{e^), (28) 

where the first two coefficients, needed below, read 

Pi =4i;(6i;2_2)-i/2, (29) 

P2=2{2E^ - E^+ 1){?>E'^ (30) 
Combined with a/M = 1 — e^, this gives 


Mn=\-h{E)e + 0{e^), 


(31) 


where b was defined in (231. Plugging (31) into Eq. (27) 
yields, in turn. 


z(o)(^.g) = 


+ 0{e^). 


V6E2 _ 2 

Finally, using (24) with (26) and (|32|), we obtain 


(32) 


^Lcons(^) ^ _ 2 )-i/2 [eH^{E-, e)] . (33) 


where gap = gap{x^', M, J) is the Kerr background met¬ 
ric, Om is taken with fixed J := Ma, dj is taken with 
fixed M, and both derivatives are taken with fixed Boyer- 
Lindquist coordinates x^. Our particular regularization 
procedure (see below) does not require the completion 
piece for r < R. The quantity is given by 


ffR — 
^compl ~ “ 


ASM) 

ffi /3 


C) 


(36) 


where the perturbations are evaluated in the sided limit 
r —> R'^ (with 0 = 7 r/ 2 ). 

Let us denote by ^Lc°recons ^''^d the contri¬ 
butions to from and respectively, 

and proceed to obtain analytically. First use 

Eq. (36) with Eq. (35) and with = g°‘^up, where 
Up = {—E, 0, 0, L} (in Boyer-Lindquist coordinates). 
This gives in terms of the circular-orbit radius 

R, its energy E and its angular momentum L. Then 
substitute the fixed-E expansions (28) for R, and [T^ 


L/M = 2E+ (6^2 - 2)i/2g 


(37) 


for L, along with a = 1 — e^. Finally, expand the result 
in e at fixed E. The outcome is 


Note that, for to be finite and generally nonzero 

(as expected) requires to blow up like 1/e for e —>■ 0. 

Our calculation of H^{E-, e) will be based on the strat¬ 
egy and numerical codes developed in Refs. [niEolEI]. 
H^{E;e) is expressed as a sum of two contribution: 


= Hi 


■^compl' 


(34) 


The “reconstructed” part Lf^cons is obtained numerically, 
starting from frequency-domain solutions of Teukolsky’s 
equation with a circular-geodesic source and retarded 
boundary conditions, following through a reconstruction 
of the multipole modes of the metric perturbation (in a 
suitable “half-string” radiation gauge [H]), and finally 
applying a suitable form of mode-sum regularization [23] 
to extract the R part of the perturbation. (A more de¬ 
tailed description will be given in the next subsection.) 
In our case of a circular-orbit so urce , the double contrac¬ 
tion of with [recall Eq. (26)] automatically picks 
out the time-symmetric part of h^p, as desired. 


The second contribution to in Eq. (34) is the 
“completion” piece which (by definition) arises 

from any part of the metric perturbation that is not 
captured by the reconstruction procedure. In our prob¬ 
lem, this piece corresponds simply to mass and angular- 
momentum perturbations of the background Kerr geom¬ 
etry (plus pure-gauge perturbations) |20l [24] . In the vac¬ 
uum region r > R outside the particle’s orbit, these sta¬ 
tionary perturbations can be written analytically, in a 
“Boyer-Lindquist” gauge, as pi] 




^0.(3 


dM 








(38) 


Notice this is an 0(e ^) quantity, so, recalling Eq. (33), 
it gives a finite contribution to We find 


(5L 


c.compl — 9 (1 ~ 


(39) 


where 


._ „-lri cons 
mpl ' I ^‘”c,compl' 

Remarkably, it follows that the completion contribu¬ 
tion, on its own, precisely saturates the censorship condi¬ 
tion (15). In the next subsection we will demonstrate nu¬ 
merically that the reconstructed part, has 

a finite (non-divergent) fixed-E limit e —>■ 0. This will 
imply 


Cj cons 
^ ^c,recons 


= 0 , 


and therefore 


ar^ = -(i-£;2). 


(40) 


(41) 


The censorship condition (15) is precisely saturated. This 
result and its implications will be discussed in Sec. [V] 


B. Numerical input 


dJ 


To validate Eq. (40), we will demonstrate numerically 
that the limit 


lini^_(E;e) =: H^{E) 

e —)-0 


(35) 


(42) 
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(taken with fixed E) exists and yields a finite value. It 
may be possible to establish this mathematically through 
analysis of the reconstructed solutions in the near¬ 
extremal near-horizon approximations (perhaps mod¬ 
elled upon the method of Ref. [IS]). Here we content 
ourselves with a numerical calculation, which, we find, 
already illustrates the finiteness of H^{E) rather con¬ 
vincingly. 

We have performed two independent numerical calcu¬ 
lations, using two different (albeit related) methods, to 
be described below. One of the methods performs best 
at e values that are not too small, and the other does 
best for e values that are not too large. The combination 
of the two methods thus allowed us access to a range of 
e values wide enough to enable taking the limit e —>■ 0 
accurately. The agreement we found between the two 
sets of results in an overlapping domain of intermediate 
e values provides reassurance. 

Our first calculation is based on method and code de¬ 
veloped by one of us (AGS, with collaborators) in Refs, 
[iniiss], with input from Ref. m- In this method, we 
first numerically integrate the sourced Sasaki-Nakamura 
equation in the frequency domain, with “retarded” 
boundary conditions at infinity and on the event hori¬ 
zon, to obtain the modes of the Weyl scalar tpQ associated 
with the metric perturbation produced by the particle. 
We then derive an appropriate Hertz potential (this is 
done algebraically in terms of ■i/'o): from which the modes 
of the metric perturbation are reconstructed by apply¬ 
ing a certain second-order differential operator [37]. We 
use a version of the reconstruction procedure that yields 
the metric perturbation in a regular outgoing radiation 
gauge anywhere in the vacuum region r > R, where R 
is the Boyer-Lindquist radius of the circular orbit. Fi¬ 
nally, we apply a mode-sum regularization procedure to 
obtain iJ^cons- The mode-sum variant we are using is 
the one developed in m, with the particle limit taken 
from r —)■ i?+. (Ref. [11] derived the regularization pa¬ 
rameter values suitable for this one-sided version of the 
mode-sum method.) The code is implemented in C++ and 
uses double-precision arithmetic. 

This has been a first implementation of the code in 
the near-extremal regime, e ^ 1. Certain technical sub¬ 
tleties arise in this regime, as recently reviewed in Sec. V 
of [25]. We have found that such subtleties were rather 
easily controlled for e > 10“^. We needed only ensure 
that inner boundary conditions were placed sufficiently 
close to the horizon and determined to sufficient accu¬ 
racy. We have not attempted to improve the performance 
of the code at lower values of e (e.g., using the techniques 
described in Ref. [3S]), but instead resorted to our sec¬ 
ond method, to be described next, whose performance 
actually improves near extremality. 

Our second calculation is based on a code developed 
by one of us (MvdM) in [21], which follows an approach 
by Fujita [28], itself based on the semi-analytical for¬ 
malism of Mano, Suzuki, and Takasugi (MST) [3^ 150] . 
In this approach, the Weyl scalar—in our particu¬ 


lar implementation—is obtained semi-analytically rather 
than numerically. “Semi” here refers to an element of the 
calculation in which a certain continued-fraction equa¬ 
tion is solved numerically. The reconstruction and mode- 
sum procedures are essentially as in the first method, but 
they are implemented using an independent code. The 
entire calculation is performed using Mathematica with 
arbitrary-precision arithmetic. 

In the MST-based calculation, working near extremal¬ 
ity is computationally advantageous. This is due to the 
improved convergence properties of the MST formalism 
for cicular orbits with a ~ I and H ~ 1/2, highlighted in 
m- In this domain, the series of special functions featur¬ 
ing in MST’s solutions for -04 converges faster. Further¬ 
more, the aforementioned continued-fraction equation is 
both faster convergent and more easily solvable (using the 
analytically-known extremal solution as an initial guess). 
As a result, the method is particularly efficient for study¬ 
ing the e —>■ 0 limit. For our purpose, it was sufficient to 
apply it in the range I0“® < e < 10“^. Above e ~ 10“"*, 
the analytically-known extremal solution no longer pro¬ 
vides an accurate enough initial guess to guarantee find¬ 
ing the solution of the continued-fraction equation for all 
frequency modes in the spectrum, making our implemen¬ 
tation of the MST formalism unreliable. 

Our calculation of H^{E) proceeded as follows. We 
considered a dense sample of E values in the range 
^'isco < E < 2. For each value of E in the sample 
we obtained a dataset i?recons > where e is sam¬ 
pled (roughly) uniformly in loge between e = 10“^ and 
e = 10“®. We switched from our fully numerical method 
to our MST-based method at around e = 10“^. H^{E) 
was then obtained via extrapolation of each of the fixed- 
E datasets to e = 0. 

For each pair {E, e} in our sample, we directly com¬ 
puted the first 70 multipoles (Z-modes) of the metric per¬ 
turbation, for use as input in the mode-sum formula. The 
remaining large-Z tail of modes was approximated by fit¬ 
ting an inverse-power-law model, as detailed in [20] . At 
high values of E, the ^-mode distribution becomes skewed 
towards larger I values, due to what may be interpreted 
as a beaming effect. A similar behavior near the light¬ 
ring of a Schwarzschild black hole was discussed by Akcay 
et al. EB, who pointed out that the implementation of 
the mode-sum technique can become problematic in that 
case, because the standard inverse-power-law tail may 
fail to manifest itself until I values larger than one can 
feasibly calculate. This effect restricted our calculation 
here to E values that are not too large—in practice, to 
A < 2. However, that should suffice for our purpose here, 
which is simply to determine the e-scaling of Hrecons 
the limit e —)■ 0 at fixed E: it is perfectly reasonable to 
assume that the e-scaling at any fixed E > 2 would be 
the same as it is for lower E. 

Figure shows e) as a function of e for 

a few E values within our sample. It is evident that 
-f^racons(-®i e) approaches a finite limit as e —)■ 0. Figurej^ 
displays the extrapolated values as a function of E. 
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We remind that the details of the function H^{E) are 
unimportant to us; we needed only establish here that 
is finite for any finite E. 



FIG. 1. Jfi^cons s-s a function of e, for a sample of E val¬ 
ues. Data points for t > 10“^ (diamonds) are from our fully 
numerical computation, while points for e < 10“'* (squares) 
are from our semi-analytical, MST-based method (there is 
an overlapping data point at e = 10“*). Error bars are 
in all cases too small to be resolved in this figure. Curves 
(dotted line) are quartic polynomial fits. At each fixed E, 
H^cons{E; e) appears to approach a constant value in the ex¬ 
tremal limit e —>■ 0. 







E 


FIG. 2. The function H^{E), obtained by extrapolating our 
numerical data for iFr^cons to e —> 0 at each fixed E. The 
actual value of H^{E) is not needed in our analysis, only the 
fact that it is finite for each finite E. 


IV. EFFECT OF STRONG FINE-TUNING 


condition would imply that overspinning is possible via 
strong fine-tuning. If, on the other hand, we can show 
that (1^ applies for any {Ei,Ef}, then censorship holds 
even allowing for strong fine-tuning. 

The evaluation of (20) requires the angular-momentum 
shift (5L);™®(i?i) and the flux ratio TZ{E). For the former 
we use our result (41). For the latter we will perform a 
numerical calculation, to be presented in subsection |IV B 


below. However, much of what we need to know about 
TZ can be deduced from simple analytic considerations, 
to be presented first. We will show that it is sufficient 
to require that TZ is bounded from below by —1/3 over 
the range Eisco ^ E < ^ m order to guarantee that 
the censorship condition (20) always holds. Our actual 
calculation will later show that TZ lies very comfortably 
above that bound. 


A. Analytical considerations 


1. Superradiant domain 


First, we consider the sign of TZ{E). We recall that 
this function, defined for circular equatorial geodesics, is 
the ratio of energy flux through the horizon to the energy 
flux at infinity. A gravitational-wave mode of the form 
~ jg Jfjjown to be superradiant if and only if 

0 < w < mflH, where fin = a/{2Mr+) is the horizon’s 
angular velocity, with r_|_ = M -|- (M^ — a^)*/^. For cir¬ 
cular equatorial orbits the gravitational-wave spectrum 
is simple: w = mfl, where H is the orbital angular ve¬ 
locity. Thus, all modes of the radiation are superradiant 
for H < Hh, giving 7?. < 0 in that case. For fl > Hjy all 
modes are instead non-superradiant, giving 7^ > 0. 

Let us now specialize to a near-extremal geometry, and 
reexpress the above in terms of a condition on the specific 
energy E of the circular geodesic. For e ^ 1 we find 
Mflir = I ~ Combining this with Eq. 


(31) translates the superradiance condition H < to 


E{6E'^ — 2) > \/2/3 (at leading order in e), leading 

to 




(43) 


as the superradiant domain in the extremal limit. Thus 

7^ < 0 for Eisco < E < Esr, (44) 

7^ > 0 for E > Est- (45) 


This will be confirmed numerically in subsection IV B[ 


We have shown that, within our first-order GSF ap¬ 
proximation, any weakly fine-tuned capture produces 
a precisely extremal geometry. Can strong fine-tuning 
push the system beyond extremality? To answer that 
question we need to evaluate the condition (20). Any 


choice of {Ei,Ef} (with Eisco < Ef < Ei) violating that 


2. Sufficient lower bound for TZ{E) 


We now show that the censorship condition (20) is sat¬ 
isfied for all Ei > Ef > ^ ii TZ{E) is bounded from 

below by —1/3. The condition involves (j) and ip, given in 
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Eqs. (18 1 and (191, respectively, where in the latter we 


now substitute for from Eq. (411. We do not know 


the sign of (j) (for given Ei,Ef) a priori, so we proceed by 
considering the two options (j) > 0 and (j) < 0 in turn. 
Case (j) > 0:— The censorship condition (20 1 becomes 


^ = + (E, - > 0. (46) 

This is trivially satisfied for Ei > 1, so it remains to 
consider < 1 , in which case the condition becomes 


<E^- x/(l - Ef)/2 =: ME^). (47) 


Recalling (21), we may bound from above using 


£+ , < 
(qc) — 


dE 


< 




s/3 


l+n{E) - l+TZn 


=: ME^), (48) 


where in the first inequality we used the positivity of 
the integrand together with Ef > and in the second 
inequality we assumed TZ is bounded from below by some 
number TZm{> —!)• Since i/i = V 2 {= 0) at Ei = 


establishing the inequality in (471 requires only showing 


that v'i{Ei) > v^^iEi) = (1 + ^ for all ^ < 1 . 

But the minimal value of over this domain is 3/2, so 
the condition becomes (1 + TZm)~^ < or 77.,„ > — |- 
We have thereby shown that the censorship condition 
(46) holds for any Ei > Ej > ^ with ^{Ei,Ef) > 0, 
under the sole assumption 


niE)>--. 


(49) 


Case (j) < 0:— The censorship condition (20) becomes 
Ip > (p'^/A, or, explicitly, 

0 < - jW;., (WJ., - 2i,{E.)) + fj., - 2E.) 
=-.A{E,,Ef), (50) 


where (p{Ei) = —( 6 iJ/—2)^/^. Since = 0 = for 
Ei = Ef, we have A{E, E) = 0 for all E > Thus, to 
establish A > 0 it suffices to show dA{Ei,Ef)/dEi > 0 
for dl\ Ei> Ef > ^. 


With the aid of Eqs. (21) and (22), we find 


dA 

[l+7^(^;,)]—=£;,+77(E,) 


^{Ei) 


- ^ 
(qc) 


(51) 

Consider the cases TZ{Ei) < 0 and TZ{Ei) > 0 separately. 
For TZ{Ei) < 0, we use ^ 4>i.Ei) (following from 

(/ < 0 ) to bound the right-hand side of (51) from below by 


i?i[l+ 377(i?i)] —2TZ{Ei)£'^^^s^. Since the last term here is 
non-negative, it is sufficient to require TZ{Ei) > — | in or¬ 
der to guarantee dA/dEi > 0 and hence A{Ei, Ef) > 0. 


If, instead, TZ{Ei) > 0, one can first use — 

[which follows from Eqs. (21) and (22), noting b{E) > 2], 
then again VV(|jc) ^ 4‘{Ei), to bound the right-hand side 
of (51) from below by Ei + TZ{Ei) [(p{Ei) + 3Ei]. This is 
non-negative for all 7?^ > ^ if and only if TZ{Ei) > — |- 
Thus, the condition ( [4^ always implies A > 0 and, in 
turn, that the censorship condition ( 20 ) holds also for 
^< 0 . 

We conclude that it is sufficient to show that the flux 
ratio TZ is bounded from below by — | in order to guaran¬ 
tee that the censorship condition ( 20 ) is always satisfied. 
In fact, recalling (45), we see that it is sufficient to obtain 
such a bound for TZ on the restricted superradiant do¬ 
main iiiisco < E < Est;. Our numerical calculation, to be 
presented below, shows that TZ is comfortably bounded 
above the value of —1/3. 


B. Numerical input 


To compute the flux ratio TZ we used our MST-based 
method described above, as implemented in [21j . The 
gravitational-wave energy fluxes to inhnity and down the 
horizon are obtained directly from the semi-analytical 
solutions for 7 / 4 , with no need to reconstruct the met¬ 
ric perturbation. Thanks to the improved convergence 
properties (already mentioned above) of the MST for¬ 
malism at e <C 1 and ~ 1 / 2 , we can obtain the en¬ 
ergy fluxes to essentially any accuracy we desire using 
arbitrary-precision computer algebra. 

To determine TZ for a given E, we calculated the ratio 
between the flux down the horizon and the flux to in¬ 
finity for a sequence of fixed-i? orbits with e values that 
decrease to 10“® in exponential steps. The value of 77(7?) 
was then found by extrapolating to e = 0. The results 
are presented in Fig. As expected, 77.(7?) is negative 
only in the range 7?isco < 7? < Egr- The minimum of 
TZ(£) appears to be attained at 7?isco with a value of 
—0.13744 ± 3 • 10“®. This is comfortably above the value 
of —1/3 required to assure that the censorship condition 
( 20 ) is satisfied. 


V. CONCLUSIONS AND DISCUSSION 

We studied the scenario in which a mass particle is 
thrown into a nearly-extremal Kerr black hole on an 
equatorial trajectory, working consistently in the first- 
order self-force approximation, i.e. taking into account 
all finite -77 effects (radiative and other) to one order in 
Tj beyond the geodesic approximation. To describe the 
fate of the post-capture geometry we followed a strategy 
set out by CB in [T^], according to which it suffices to 
consider two types of captured orbits near the capture- 
scatter separatrix: (i) weakly fine-tuned orbits, which ex¬ 
ecute 0 (ln? 7 ) quasi-circular revolutions below the ISCO 
prior to falling into the black hole, radiating 0(r/^ In 77 ) of 
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FIG. 3. Numerical values for the flux ratio TZ, as a function of 
specific energy E, for unstable circular equatorial geodesic or¬ 
bits in the extremal Kerr limit. Each dot represents a numer¬ 
ical measurement, with error bars being much smaller than 
the size of the dots. Orbits with E < Ebt = 2/-\/3 are su¬ 
perradiant, with TZ < 0. The inset expands the area around 
the minimum at Eiaco = 1/y/Z. We find a minimum value of 
TZ = —0.13744±3-10~®, safely above what is required for the 
censorship condition (|20[) to hold (dashed red line). 


gravitational-wave energy in the process; and (ii) strongly 
fine-tuned orbits, which execute 0{ri~^) revolutions and 
emit 0(r?) of energy. 

In Sec. in we found that, within our first-order GSF 
approximation, any weakly fine-tuned capture leads to 
a precisely extremal geometry. This implies [12] that 
“generic” captures (ones that are not fine-tuned at all, 
with L — 2E negative and not small) produce sM&-extremal 
geometries. In Sec. IV we further established that strong 
fine-tuning promotes censorship: all strongly fine-tuned 
captures produce subextremal geometries. Thus, one can 
at best reach extremality, using weakly fine-tuned orbits 
(any such orbit would do), but there is no way of over¬ 
spinning the black hole. In summary: 


Within the first-order self-force approxima¬ 
tion (and excluding deeply bound orbits), 
equatorial captures generically result in a 
subextremal post-capture geometry. One can 
at best achieve extremality, through weak fine- 
tuning, but overspinning is not possible. 


That overspinning appears to be possible in the geodesic 
approximation [T] is simply an artefact of ignoring im¬ 
portant GSF terms that appear already at leading order 
in the relevant overspinning conditions. 

We note that the above conclusions were arrived at al¬ 
most entirely via analytical considerations. We required 
only two pieces of numerical input, one confirming the 
boundedness of the extremal limit in Eq. (42), and an¬ 
other establishing the bound (49) for the flux ratio. Both 
numerical computations involve only circular geodesic or¬ 
bits, and neither requires particularly high precision. The 
above general conclusions are entirely robust with respect 


to numerical error. 

However, it is important to remember that here we 
have been working strictly within the framework of the 
first-order GSF approximation, with no control whatso¬ 
ever over high-order GSF corrections. Since the first- 
order analysis appeared to allow for an exact saturation 
of the overspinning condition, higher-order effects may 
qualitatively change the outcome. A second-order anal¬ 
ysis may potentially yield any possible result: that the 
final geometry is always subextremal, or that overspin¬ 
ning is possible, or (once again) that the black hole can 
at most be brought to extremality. In that respect, our 
first-order GSF analysis—just like the geodesic analysis 
of Ref. 0/ —does not provide a conclusive answer to the 
question of overspinning. It is not clear if the question 
can be fully resolved at second-order or at any other finite 
order in perturbation theory. This may be a disappoint¬ 
ing conclusion, but it is an interesting one nonetheless. 

For whatever it is worth, let us return to discuss the 
consequences of our first-order analysis. We have found 
that overspinning is not possible, consistent with the con¬ 
jecture of weak cosmic censorship. However, we have also 
founds that, through (weak) hne-tuning of the orbital 
parameters, one can drive the black hole to extremality. 
This is an intriguing possibility, because it appears to 
be in violation of Israel’s “third law” of black hole dy¬ 
namics [32]. We have not studied in detail whether the 
conditions of the third law can be said to be met in full 
in our problem. Reference [7| contains some discussion 
of this point, and proposes how any apparent violation 
of the third law in the overspinning problem might be 
reconciled. 

There are several ways in which our analysis may be 
improved and further tested. First, it would be desir¬ 
able to repeat the calculation of the angular momentum 
shift using a direct integration of the GSF, via 

Eq. ^ (with Fa —>■ E“™®). This would eliminate our 
reliance [in deriving Eq. (24)] on the effective Hamilto¬ 
nian formulation of Refs. HlHISI, which is axiomatic in 
nature. In fact, an explicit demonstration of agreement 
between the direct formula (j^ and redshift formula (24) 
would constitute an important test of the Hamiltonian 
formulation. There is work in progress at Southampton 
to directly (numerically) evaluate via Eq. (j^. 

Second, it may be possible to replace some (or all) of 
the numerical input for our analysis with analytical argu¬ 
ments. Specifically, one may se ek to establish analytically 
the boundedness of in Eq. (42), and the lower bound 
(49) for the flux ratio TZ. Both may be achieved by ex¬ 
tending the ISGO analysis of Ref. to unstable circular 
orbits, and (in the case of H^) also from the Weyl scalar 
to the reconstructed metric perturbation. We leave this 
to future work. 

To conclude, we reiterate our view that our work 
represents a first complete analysis of the overspin¬ 
ning/overcharging problem through first post-geodesic 
order in perturbation theory, for a particular capture sce¬ 
nario. Of course, we have only explored here a fraction of 
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the space of interesting scenarios. Our analysis did not 
cover, for example, (i) very low energy configurations of 
deeply-bound orbits, (ii) non-equatorial orbits, (hi) ul- 
trarelativistic or null particles, or (iv) spinning and/or 
electrically charged particles on a Kerr-Newman back¬ 
ground. These scenarios all deserve examination. 
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